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Abstract 


Self-consistent dynamics of the multi-component, A^-field quintessence and 
gravity is formulated as relativistic A^-body dynamics, embedded in a static 
viscous flat space and under the forces given by the interacting Lorentz scalar 
potential via exchanges of field bosons. The Ratra-Peebles power-law po¬ 
tential of effective single-field quintessence can be derived from this ’’micro¬ 
scopic” perspective. In certain situations, the effective dynamics can be made 
identical to that of single complex quintessence, except for that the overall 
U{1) symmetry is not manifestly broken. The present formulation provides 
a convenient gauge for analyzing superhorizon perturbations and possibly for 

quantization of superhorizon fields and gravity together. 
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I. INTRODUCTION 


The discovery of accelerating Hubble expansion has inspired intense investigations on 
various acceleration mechanisms [^]-0]. Conventionally, the cosmological constant of a suit¬ 
able value does just fine to explain satisfactorily all existing data However, as the 

quality of cosmological observations are rapidly to be further improved in the near future, 
it is therefore timely to explore possibilities beyond what has conventionally been taken for 
granted. The proposal of quintessence marks an important breakthrough in how one per¬ 
ceives the major constituents of the present-day universe, which have conventionally been 
believed to be describable by classical physics. This new constituent is a field that permeates 
the entire space in a relative uniform manner, and is not a collection of matter particles. 
Investigations on the quintessence dynamics and its ultimate cosmological implications have 
been an active research area in recent years |0]. 

The field description in the context of cosmology is not new, but was in the past limited 
to the physics of very early universe, where classical physics no longer holds. The inflaton 
field that drives the early-universe inflation is a pronounced example. Reheating is another 
example; it addresses how the radiation is created by dissipating the rapid oscillation of the 
same inflaton field previously driving the inflation In fact, the proposal of quintessence 
is well-motivated, by the recognition that the origin of the cosmological constant can be the 
ground state of a certain scalar field. If so, it is not too far-fetched to consider the field to 
be now evolving on its way to settle to the ground state on the cosmological time scale. 

Somewhat different from the pure scalar field dynamics in the very early universe, dy¬ 
namics of quintessence must generally take into account other evolving energy components, 
such as matter and radiation, in the universe. The presence of other energy components 
turns out to dilute the strengths of both the ’’drag” resulting from Hubble expansion and the 
’’force” given by the quintessence potential. It is therefore of relevance to investigate how 
all coupled energy components, including gravity, evolve in time under a unified framework. 
This paper first aims to address this issue. 
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On the other hand, there may be more than one quintessence helds present in the uni¬ 
verse. In general, these N helds couple through some mutual interactions, in addition to 
the gravity. The complexity of the multi-held problem rises drastically, due partly to the 
increasing number of N degrees of freedom, and also to the N‘^ mutual interactions. Nev¬ 
ertheless, when the number of helds N is sufficiently large the problem can become simple 
again, since one may extract a few relevant dynamical variables and derives an ehective the¬ 
ory from this many-held problem. In an analogy to many-body problems, such an ehective 
theory is derived from a held theoretical approach. If this program can be made successful, it 
will provide good foundation for the ehective potential of the ehective quintessence in terms 
of elemental mutual interactions of individual quintessence helds. Furthermore, the mutual 
interactions should be regarded as exchanges of some mediate bosons in the held space. The 
present paper mainly aims at formulating a useful framework for such a many-held problem. 

It also turns out the present formulation, using T(= ln(a)) as the new time variable, is 


most convenient for analyzing cosmological perturbations in the long-wavelength limit [10 


In this gauge, T can remain unperturbed to the order of k/H, the ratio of Hubble radius 
to the wavelength of superhorizon perturbation, and the cosmological perturbations can be 
conveniently perceived as the evolving diherence of two homogeneous and isotropic universes 
with slightly diherent initial conditions. 

In this paper, we do not wish to distinguish the quintessence from the scalar helds that 
were active in the early universe before photons were produced (reheating), and regard these 
helds al together as the quintessence. That is, one may turn oh the photon huid in the present 
formulation to address the inhaton dynamics. The dynamics of single quintessence held is 
formulated in Sec.(2). Extension to multi-component, multi-helds is presented in Sec.(3), 
where examples for deriving the ehective quintessence models are given. We also touch 
upon the cosmological perturbations with the present formulation in Sec. (4). Conclusions 
and comments are given in Sec. (5). 
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II. SINGLE-FIELD QUINTESSENCE 


Assuming a homogeneous and isotropic universe, we have the equation of motion for the 
scalar held 0 under a potential V (0) as 


0 + 3hf0 -|- 


0U(0) 

(90 


0 . 


( 1 ) 


The Einstein equations reads 

where g = AnG/S, 0' = d(()/dT with T = Ina, 10.// = U(0) + pf, and pf and pf are the 
energy density and pressure of other huid components, which may also include the space 
curvature. The space curvature can generally be regarded as a huid in this homogeneous 
and isotropic limit, with an equation of state Pf = —pf/3. The time dependence of p/ is 
proportional to and its value can be either positive or negative depending on whether 
the space has negative or positive curvature, respectively. 

These coupled equations can be simplihed if we let the independent variable of Eq.(l) 
be T. Thus, 0 = H(()' and 0 = if^0" + Hep', with H and H replaced by Eq.(2). Equations 
(1) and (2) combine to yield a single equation: 




(3) 


2K// ' 2g .90 

The peculiar feature of this equation is that y^0' = ±1 are the hxed points, and solutions 
with p0'^ — 1 < 0 or > 0 cannot cross the hxed points to change the signs. The peculiar 
behavior at the hxed points means that there exists a maximum rate of coherent change of 
0 over the horizon. This condition is also rehected in the expression of in Eq.(2), where 
may become singular. If 10// does not change sign, the singularities in can not be 
removed by having the numerator to vanish simultaneously with the denominator. Such an 
unphysical result is automatically avoided in Eq.(3), demanding that p0'^ — 1 must remain 
of the same sign. As a non-negative 10// is of physical relevance, we shall focus on the causal 
regime where p0'^ < 1 in the rest of this paper. 
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Having understood so, we may make an analogy of the field dynamics to particle dynam¬ 
ics, in that one identihes ^/g(|)' as the velocity w, (1 — as the Lorentz factor 7 , and 

u = as the momentum. This nonlinear transformation puts Eq.(3) into 


du 

dT 


+ 3(1 


PL±PI)u + ^ >9 In He// 
2Veff ^ 2^ d(t> 


0 . 


(4) 


The hrst and third terms are the ”01” components a second-rank tensor and the second term 
is a vector. It demonstrates that this dynamical system is not Lorentz covariant, because the 
frictional force given by the Hubble expansion behaves as if the system contains a viscous 
medium which has already chosen a preferred reference frame. 

Of particular interesting of Eq.(3) or (4) is the softening effect to the held 0 provided by 
another energy component p/. The held has a vacuum state 0^, where H(0„) = = 0. 

The ’’force” —(91n Veff/d(j) can be locally expanded as —2(0—0mm)[(0—0mm)^ + (2p//H")]“h 
It is analogous to a two-dimensional Coulomb attractive force given by a hnite-sized charged 
distribution. If the soft core does not exist, i.e., p/ = 0, the ’’kinetic energy”, or 7 , will 
become inhnity at the vacuum, so that in Eq.(2) remains hnite. 

To put the above into a quantitative perspective, we note that rs-/ Epianck, the Planck 
energy. Hence the dimensionless ’’displacement” ^/gcf) and velocity ^/gcf)' are all normalized 
to the Planck scale. The Planck scale thus typically implies a ’’relativistic” regime. However, 
the reverse statement is typically not true. Energy scale much below the Planck regime can 
also be ’’relativistic”, as long as the dynamical time scale is sufficiently small so that d/dT 
is large despite 0 being small. Such a situation is commonly encountered when Veff 0. 

Despite not Lorentz covariant, the dynamics can still be put in a Lagrangian framework 
that describes the self-consistent evolution of held and gravity. A straightforward inspection 
shows that the action for the coupled dynamics becomes 


I LdT = -I e^^dTr^pgVeff (5) 

By dehning a new potential <h(0,T) = InVe/f + 6 T, we have here an action of particle 
dynamics in the presence of a dilation held: ^/2^ f exp(<h(0, T)/2)(ir, where <I> is a dilaton 
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field and r is the proper time. In fact, the action given in Eq.(5) can be derived rigorously 
by summing up the Lagrangians of quintessence held, gravity and other energy components, 
as will be described in the next section. 

Although there is no obvious conservation law for 0 due to cosmic expansion, this La- 
grangian nevertheless allows us to obtain a Hamiltonian, which can yield an approximately 
conserved quantity. This Hamiltonian h is constructed by the usual procedure 

h = <P\dLld<P') -L = ^^2gVeffe^^. ( 6 ) 


Actually, h = ifexp(3T) according to Eq.(2). 

To show h to be indeed useful, we obtain the equation for h by multiplying Eq.(5) by u 
to get 


d\nh 3 pf + pf 
2Veff 

or after a straightforward algebra 


d 

dT 


{V + 


-Qu^Veff. 


(7) 


( 8 ) 


It is clear that Eq.(7) is appropriate for the ’’ultra-relativistic” regime, where the right-hand 
side can be ignored and one obtains approximate conserved quantities. 

The ’’ultra-relativistic speed” means that ^/g4>' —±1, implying ^/gcf) ±T in the large 
|T| limit. The ’’Lorentz” factor 7 can stays large as long as 2gVeff << /i^exp(— 6 T). On the 
other hand, since h = const, the Hubble parameter H{= dT/dt) is proportional to exp(—3T) 
in this regime. Solving for T{t){= dln{a)/dint), we hnd that exp(T) ~ a{t) oc t^^^, i.e., an 
universal expansion regardless of the form of V ( 0 ) in this kinetic-energy-dominated regime. 
This is the slowest asymptotic Hubble expansion rate one may hope to obtain. 

From Eq.(4), one also sees that apart from the usual notion of inflation where a large 
range of flat V{(j)) and a potential minimum are needed for it to occur, the inflation can 
also take place in a peculiar situation when pf = 0, the potential has a maximum, and the 
incoming kinetic energy is tuned to have just the right value for it to be completely consumed 
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upon climbing to the potential top. Since u = 0 and (91nl//(90 = 0, a state situated exactly 
at the separatrix also gives a constant H, thereby yielding inflation. 

Before leaving this simplest single-field dynamics, we shall also look into the dynamics 
of typical relevant potentials. 


A. Exponetial potential 

The exponential potential, V{(j)) = exp(6/5y^0), has interesting conformal symme¬ 
try, when pf = 0 and (3 = ±1. The conformal transformation, = (1 ± u^/l =F 

- VtT) - To and f = (1 ± Vt/1 =F - Vt^cj)) + Tq - (1/6) ln(l ± 

Vt/1 =F Vt), transforms the original action / dT^l — 5f((i0/(iT)2exp[3(±y^0 -|- T)] into 
/ dT^l — exp[3(±y^0 _l_ where the transform ’’velocity” Vt and the shift 

To are both constants. The transformed action assumes the same form as the original ac¬ 
tion, and thus yields the same equation of motion. In fact, this particular case |/3| = 1 lies 
at the boundary of two distinctly different regimes, to be discussed below. 

When l/5| > 1 , it is possible to make the above transformation such that the T de¬ 
pendence vanishes by choosing Vt = —This is a space-like regime where the La- 
grangian becomes time-independent and free of friction. The conserved Hamiltonian h 
is a constant of motion and becomes h oc 7 exp( 3 y^ 0(/5 ± 1)). It allows for a solution 
0 = [l/3^{f3 ± 1)] ln[cosh(3(0 ± l)(f - n)/2)]. 

On the other hand, when \j3\ < 1, one may choose Vt = —j3 so that the original exponen¬ 
tial potential is transformed to become a 0 -independent vacuum energy, thereby yielding 
solutions that experience only the frictional force. This is a time-like regime. The La- 
grangian becomes oc 7 “^exp( 3 (l ± j3)T). Since there is no applied force in the transformed 
coordinates, we have an asymptotic solution at rest, 0 = 0 and dcjy/dT = 0. It then follows 
that H oc 7 exp( 3 / 5 y^ 0 ) oc exp(3/5utT) = exp(—3/5^T). Hence, the scaling factor a oc 
For 0^ < 1/3, these solutions result in power-law inflation, for which the universe undergoes 
accelerating expansion. 
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B. Power-law potential 


The repulsive power-law potentials, V{(j)) = for n > 0, can also yield inflation¬ 

like solutions, known as the Ratra-Peebles potentials [n]. The inertia, i.e., the second 
derivative of Eq.(3) or (4), can only be relevant only for a short while for the field to pick 
up the ’’speed”. But it will quickly become subdominant, and the field reaches the terminal 
’’velocity”, where the frictional force balances the applied force. Moreover, as the terminal 
speed is always ”non-relativistic”, Eq.(3) or (4) gives the solution, cj) oc and hence 
H oc y/V oc T or a oc exp(t^d2+n)^^ inflationary solution. The above is for the 

case pf = 0. When p/ is finite and even dominates K//, fhe balance of friction and applied 
force then gives E(0) oc where we have let pf oc a~^. Apparently V declines more 

slowly than pf does, and hence the scalar field will eventually dominate and the Hubble 
expansion approaches that of the fluid-free case. 

The dynamics can be more complex when the power-law potential is attractive, V (0) oc 
0^"^, as it exhibits anharmonic damped oscillations. Nevertheless, the following general prop¬ 
erty holds. First, since the attractive force, —01nl4///50, becomes increasingly stronger 
in the regime where the average potential strength {V) » pf, the oscillation frequency 
is increasingly higher. When the oscillation frequency is sufficiently high, one has a small 
number, the ratio of the oscillation period to the Hubble time, as an expansion parameter 
for a perturbative treatment in evaluating how the oscillation amplitude declines in time. 
(This high-frequency regime may also occur even when pf » {V) because pf can decrease 
rapidly in time. An accurate assessment of when the high-frequency regime holds will be 
deferred to the end of this subsection.) 

To the leading order, the orbit in one oscillation can be viewed as dissipationless, thereby 
allowing one to calculate the ’’ideal” orbit. One may subsequently substitute the ’’ideal” 
orbit into the damping term to evaluate how much damping it incurs to the orbit in one 
oscillation period. The long-term damping of the oscillation amplitude can therefore be 
evaluated to the leading order after these two steps. 
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Ignore the right-hand side of Eq.(8) and one gets the momentum of the dissipationless 


oscillation orbit: 


^2 ^ E(T) - V(,P) 


(9) 


V{<l>)+p,(TY 

where E is a quasi-static constant of motion, slowly dissipated by the right-hand-side of 
Eq.(8). When pf >< V >, the dynamics assumes non-relativistic motion and hence E ^ 
Pfv'^ + 1^(0). That is, quintenssence acquires an effective mass of 2pf. Averaged over one 
oscillation period, Eq.(8) gives, to the leading order, the rate of secular change for E: 

DE fdTiV + PfW j d(pxlE - V[(p) ^ dE 


^JdT{V + pf)u^ ^ Jd(t)^jE-V{(j)) 

-D_ __ _= —6- 


= -3- 


( 10 ) 


DT “ IdT Jd4,{pE-V{4,))-^ dMJd4>^jE-V{4,)) 

Here we have used the non-relativistic approximation dT = y/gd(j)v~^, D/DT is the deriva¬ 
tive on the slow time scale and the second equality is merely an algebraic equality. Since 
the right-hand side of the third equality is a derivative on the slow quantity, Eq.(lO) can 
thus be trivially integrated to yield 

1 


T{E) - To = -- ln(| #^T-E(0)), 


( 11 ) 


where Tq is also an integration constant. 

We now apply the above results to the attractive power-law potential, 1^(0) oc 0^”. 
The integral /— V (0) oc and T(Vo) = — lnT[(n -|- l)/6n] -|- const., or 

E cx That is, a quadratic potential, n = 1, yields E cx a~^, similar to a matter 

fluid and a quartic potential, n = 2, gives E cx a~‘^, similar to a photon fluid. 

Having Eq.(ll), we may further evaluate Dt = J dt{= f dTH~^), so that the averaged 
Hubble expansion rate H = DT/Dt has the expected expression: 

^ = ^2g(E(T) + p(T)). (12) 


Equations (11) and (12) together permit the deceleration parameter to also be evaluated: 


d{T) 


1 D^a _ i citi(g(r) + pj(r)) 

2aH^ 2 DT 


(13) 
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All the above results for the long-term evolution is based on the assumption that oscilla¬ 
tion frequency is much larger than the Hubble expansion rate. It is evident from Eqs.(3) or 
(4) that when p/ is much greater than V (0), the force strength is substantially reduced by the 
factor V (0) /pf, thus reducing the oscillation frequency by (0)/p/. Hence, there is a limit, 
beyond which the above results fail to hold. The regime in which this may happen is always 
” non-relativistic”, and one may approximate the right-hand side of Eq.(4) by —where 
00 is the dimensionless oscillation frequency, to estimate oo [{dV/d(j))/2^(j){pf + 

Again, using the power-law potential, V oc 0^”, we hnd that oo oc 0“^ oc if {V) » p, 

and that u oc qc when p » (V), where b = 3(1 + w)/2 — 3{n — l)/(n -|- 1) 

and the fluid equation of state w = Pf/ Pf has been assumed. In the former case, all conhning 
(n > 0) potentials give the desired high frequency. In the latter situation, we, however, need 
6 > 0, or n < (3 -|- tc)/(l — w). For photon fluids, one needs n < 5, and for matter fluids, 
n < 3. 


III. DYNAMICS OF MULTI-COMPONENT N FIELDS 

A. General framework 

The success in describing the field dynamics by relativistic particle dynamics motivates 
us to go one step further, and extend this formulation to the dynamics of multi-component, 
multi-helds. Here, the component of a held refers to the requirement that it has a geometric 
structure as a vector. 

As an illustration for how this can be done, we shall hrst consider a simple example of 
multi-component, single-held quintessence with no huid component present. In this case, the 
potential U is a scalar and cannot be an arbitrary function of 0*, where i is the component 
index. A particularly interesting potential is H(0) = Voexp[0 • 0], where f3 satishes 0^ = 1. 
In this particular case, the constant vector 0 picks out a particular direction of 0. Along 
this direction, one may perform the conformal Lorentz-t-shift transformation discussed in 
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Sec. (2) for the single-component quintessence, and obtain an invariant action under such 
a transformation. The action is also invariant to the conformal helical (rotation-l-shift) 
transformation along the (3 direction. 

It turns out that the action of this particular iV-component, single-held quintessence 
can also be the action of the N — M component, M-£eld quintessence. This is due to 
that the Lorentz factor (c.f., Eq.(5)) in the action contains the ’’velocities” of all helds, as 
will be elaborated below. Hence, the above conformal symmetry still holds for the N — M 
component, M-£eld quintessence, as long as (ff = 1- 

The above example shows that when the potential V assumes some particular forms, 
there can be no distinction between the multi-component single-held and the multi-held. 
To have a distinct diherence, the potential should have a more general form. Below, we 
consider such a situation. For the purpose of illustration, we are conhned to the two- 
component helds, which can be conveniently described by complex helds. Fields with more 
than two components can be easily generalized. 

Let xjji be the i-th of a system of N two-component helds. We are interested in an 
autonomous system of interacting N helds, where the interacting potential can be written 
as Vi^j = F(|'0j — V’jI) with i ^ j. The action for the system of N helds, the gravitation and 
a huid component can thus be written as 


‘S'o - J 

i ij y 


+ (A( 


n 


A 


Pf+Pf 


2{pf+pf) 


)+Pf)l (14) 


where is the space volume, R{= (1/2)[(a/a)^ -|- (a/a)]) the scalar curvature, A the La- 
grangian multiplier, and nf and 6f are the number density and velocity potential of the 
huid, respectively. 

The Lagrangian multiplier A is introduced to dehned what the huid 4-velocity is meant 


to be. In fact a full huid Lagrangian has been explored before []^, and here we consider 
only the spatially-uniform limit. Variation of So with respect to A gives the constraint: 


^Pf+Vf 


1 ; 


(15) 
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variation with respect to 6f yields a conservation law, 

_ 2 


= s 


( 16 ) 


Pf+Pf 

where s is a real constant. Since the 4-velocity satisfies = 1, the constraint, Eq.(15), 

defines the fiuid 4-velocity as = Ufd^Of/{pf +p/). It then follows that Eq.(16) gives 
a^UfUo = as the fiuid is rest in the comoving-frame of cosmic expansion, this becomes 
a^Uf = i.e., the conservation of particle number density in an uniformly expanding fiuid. 

Variation with respect to Uf gives d{\pf -|- (A — l)pf)/d\nnf = \{pf +Pf)- Note that the 
ideal-gas law demands dpf/dlnuf = pf + p/, therefore fixing A = 1. 

The Einstein equations is derived from this action by variation with respect to the scale 
factor a, and it gives, as before. 


tt2 _ _ 2 J Vi,j + Pf Pf _ fi 


-2H^ + ‘ig(2Y,V,,,+p,-p,), (17) 


where = H~^ipi and 

Now, replacing dt by dT/H, substituting H of Eq.(17) into d/dT and recognizing R = 
6g{ptotai — ^Ptotai), we finally obtain the effective iV-field and fiuid action for Eq.(14), much 
like how Eq.(5) is obtained. In this effective action the gravity is manifestly hidden: 



where (T) ^ “ gn%0'ff/{pf + p/). 

The action recovers Eqs.(15) and (16), and again the appropriate ideal-gas law is needed 
to fix A = 1. The non-trivial extension of this action from that of Eq.(5) is the appearance 
of an effective Lorentz factor f, which is not the sum of individual Lorentz factors as one 
may have expected. Rather, it is the sum of ’’kinetic energy” of all fields that contributes to 
T. When N is sufficiently large, it is obvious that << 1 in order for T to be real, thus 

demanding each individual field to be ”non-relativistic”. In addition, it is straightforward to 
derive the corresponding Hamiltonian h, as we did for Eq.(6), and the relation h = H exp{3T) 
still holds. 
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Although Eq.(18) contains a complete Lagrangian for describing the coupled dynamics of 
fluid and quintessence in terms of the logarithmic scaling factor T, the fluid and quintessence 
in fact form a master-slave system. That is, in terms of the variable T, the fluid dynamics 
is independent of the presence of quintessence, but the quintessence dynamics is governed 
by the presence of fluid. As the fluid dynamics can be easily solved when the equation of 
state is given, it is therefore possible to simplify the algebra by focusing on the quintessence 
piece of Eq.(18), with the fluid variables pf and pf treated as known functions of T. This is 
the strategy that was adopted in Sec.(2) and it gives a much simpler and practical result: 

S, = -/3 J e^^dTr-^^2gVeff, (19) 


with the fluid satisfying 


dpf 

dT 


= -3(p/ +P/), 


( 20 ) 


where V^ff = J2i,j ^j+P/(^) T”^ = — g Moreover, the Einstein equations, 

Eq.(17), now become Eq.(2), with the ’’Lorentz” factor in Eq.(2) replaced by the multi-field 
T defined here. 

We would like to mention in passing that the way by which the multi-field version of 
Eq.(3) is derived from this multi-field quintessence action, Eq.(19), has a slight twist. Extend 
Eq.(3), G = 0, to the expected multi-held version, Gi = 0. The Euler-Langrange equation 
of Eq.(19) in fact yields Ej[(<^p - Pij) + ^~^Pij]Gj = 0, where Pij = dij - Efc(V'D^ 
a projection operator perpendicular to V’)- Since {6ij — Pij) + V~‘^Pij is non-vanishing, this 
equation simply implies that the component of Gi parallel to V’) and those perpendicular to 
ip) both vanish, and therefore Gi = 0. The projection into two components relative to the 
’’velocity” manifestly indicates that the special-relativistic effect is in action for producing 
anisotropic forces. In particular, Eq.(4) should be rewritten from du/dT = F to dui/dT = 
~ Pij) + ^~‘^Pij]Pj in the multi-held case. This feature of dynamics given by Eq.(19) 
has been discussed extensively in the past, and the potential V has been known as the 
Lorentz scalar potential . 


13 






B. Coulomb field gas 


One may be contented with this Lagrangian, where Vi^ is a given function. But one may 
also be somewhat more ambitious and asks whether there may be a fundamental dynamical 
equation that accounts for fdj. Actually, such a question is well-motivated, since these 
A^-fields are like A^-particles. It is natural for A^-particles to interact via exchange of some 
mediate bosons, such as the familiar photons in electromagnetic interactions, and these 
mediate bosons obey their own dynamical equations. Here, we shall regard Idj as the 
manifestation of such scalar mediate bosons of zero frequency. To attack this problem, we 
will need to go beyond the original discrete field -ipi and consider V to be the field generated 
by and interacting with the discrete '^j, which in the continuum limit can be replaced by 
/ where the Dirac ^-function is in use. The simplest version for V(ip) dynamics 

has an action: 


where is the coupling constant with e having a dimension of energy, and k is also a 
constant having a dimension of energy as well. 

A good working hypothesis for the constant k, which relates the metrics of 'ip to the 
expansion factor T, is that k ~ thus governed by the physics of gravity. On the other 

hand, the field-field coupling constant e'^ is likely governed by physics other than gravity, 
and we shall leave its magnitude as a free parameter. 

Moreover, one needs to extend the action S'o of Eq.(14) to a more general one: 


S 


^3T 




:(^i) 


R 

l^g 


2 P/+P/ 


iPf + Pf))+Vf] + Sv, 


j 


( 22 ) 


Variation of S with respect to V gives the equation of motion for V: 


e d e^'^ dV 



(23) 
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Equation (23) shows that there exists an retarded effect in V, and hence the functional 
form V{\ipi — t/’jl), suitable only for static helds, is over-simplihed. However, as has been 
mentioned earlier, each individual held generally has a ”non-relativistic” velocity in 

the large-iV limit and moreover ~ g~^. Therefore, to the leading order of the held 

V is approximately quasi-static. This is reminiscent of the situation where a collection of 
slowly moving particles of same charges generate mostly the electrostatic huctuations and 
only little electromagnetic huctuations. Ignoring the time derivative of V in Eq.(23), we 
obtain a static, axi-symmetric solution for the potential contributed by the discrete source 

(24) 

or the Green’s function for the two-dimensional Poisson equation. 

Finally, we sum up all terms in S explicitly to obtain the total ehective action. It follows 
after an integration by part over ip that the total action is almost identical to Sq of Eq.(18), 
except for replacing Vi^ by (l/2)Id,j- This factor 1/2 is characteristic of a self-consistent 
electric or gravitational potential built up from a vacuum. The contribution of interaction 
potential to the Einstein equations is also given by the new potential, (1/2)1/^, rather than 
the original one. The above shows that quintessence, its interactions, gravity and huid can 
all be incorporated in a single action. 

But, as has been mentioned earlier, the full action with the huid included has no obvious 
advantage for algebraic manipulations, and hence we shall from now on remove the huid 
action from the quintessence action. Equation (20) will be needed for providing p{T) that 
appears in the quintessence action. Repeating the same procedure, we hnally arrive at an 
explicit form for the quintessence ehective action, Eq.(19): 

= e^'^dTr-^ 2g[^Y.^n{\ipi-ipj\) + pf], (25) 

from which the equation of motion for each individual held ipi can be derived. 

The above exercise demonstrates that the quintessence potential V can be given a more 
fundamental origin, and in this particular example we obtain a system of N helds interacting 
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among themselves like a two-dimensional Coulomb gas. If the helds have three components, 
the interacting system will be like a three-dimensional Coulomb gas. Generally, the m- 
component helds will behave like a m-dimensional Coulomb gas. The Coulomb gas self-repels 
and can maintain a fairly uniform density. It can therefore be a relatively simple N-body 
dynamical system for the description of the cosmic evolution, if the Coulomb gas is initially 
cold and uniformly distributed. For a hnite N, the uniformly (in coordinate) distributed 
helds occupy within a uniform disk of hnite area vrg^. The disk experiences a uniform radial 
force, except at the edge. The outward force at the edge, —dV/dq = — 
proportional to characteristic of a two-dimensional Coulomb force. The radius q thus 
becomes the only relevant dynamical variable for describing the evolution of the held energy 
density, and we recover the single-held dynamics. 

Specihcally, Eq.(25) contains T”^ = — g and J2i,j Vij- In a uniform cold disk, 

the velocity is proportional to the radius much like the Hubble expansion, and hence 
= {d\nq/dT)'^{27r \i/jfd\i/j\){N/7rq‘^) = {N/2){q'y. On the other hand, the force 

at the edge can be obtained easily by using the Stokes theorem for a disk of uniform surface 
density e'^/vrg^: that is, V{q) = e'^ln(l/g) -|- cq, where Cq is a constant serving to ohset 
the ground-state energy. The ehective action for the single-held description of this A^-held 
problem now becomes 

or upon being properly normalized. 


Sq,eff = - /’e^^dT(l - c/(g')^)i/2[2c/(e^ln(^) 

J u 


(27) 


where q = q^N/2 and the constant Cq is so chosen that the maximum q is unity, a require¬ 
ment that the initial condition must satisfy. This action describes the m = 2 case of the 
Ratra-Peebles Potentials |^, V{q) cx for the single-held quintessence, where late 

inhation can always occur. In fact, all Ratra-Peebles potentials of integer m — 2 can be 
generated by this A^-held model of m-component helds. 
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C. Warm Coulomb field gas 


The above cold-gas model is of the simplest type. A reasonable extension of it can be 
for an initial condition where all helds are still uniformly distributed but having thermal 
velocities. In this case, we need to consider a mean-held phase-space distribution fiip, ip') to 
describe the Wheld kinematics when N is large. That is, while V(q) remains to be —e^ ln(g) 
as is for the cold gas, the quantity demands a more elaborateed treatment. In the 

large N limit, the long-range nature of Coulomb interactions allows each held to have a 
large mean-free path and ’’sees” only the mean held Vmf, which equals to V{q) at q = q. 
According to Eq.(3), the distribution function therefore satishes 


(9T ^ ^ dip 


-T 


- 2 [ 


1 din Vmf _ _ P/+P/ ^^./j . 9f 


^2g dip Veff 

We follow the mean-held relativistic kinetic theory: 


dip' 


= 0 . 


(28) 


1 ^ = /crv>/= jd"',pni„ (29) 

i 

where m is the number of held components and the number density n?, has been normalized 
to unity; 

P E = / d'ii I i’’)] = / (so) 

i 

where Vb is the local bulk velocity; 

= J d'^'ip j (Tip'lgiip'ffiip^ip')] = J d'^ipiubvl + (31) 

i 

where is the thermal energy. Here, we have assumed that / has an isotropic velocity 
distribution. 

The appearance of the thermal energy adds an extra dynamical variable to this 
problem, and it needs to be solved self-consistently. This task can be relatively easy to 
handle when N is large. To a good approximation, each held freely streams, much like an 
ideal-gas particle, so that the adiabatic ideal-gas law provides the needed description for U^. 
For uniformly distributed A^-helds, we may use a scaling factor b{t) to account for the bulk 
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motion, much like how the scaling factor of the Hubble flow a{t). That is, ipi = b{t)xi and 
-ipi = Hip' = hxi + bxi = H{vi + 6vi)/^, where 5vi is the ’’thermal” velocity. In fact, the 
individual held tpi satishes the original held equation, 

0 = ^ = ib + 3Hb + iv:;^f/2)b]x, + [bx, + x,{2b + 3Hb)]. (32) 

U'lpi 

The hrst bracket vanishes, as it describes the bulk how or the uniform expansion. The 
second bracket must therefore also vanish, thereby allowing Xi to be solved: Xi oc b~‘^a~^. 
Furthermore, the dynamical variable, the disk radius q, is proportional to the held scaling 
factor b, and it thus follows that the thermal energy 

= (b II (33) 

i 

where a is a real constant. 

Since bvf has already been solved, there is no need to consider the ’’thermal” compo¬ 
nent of g{ip'P)‘^ appearing in of Eq.(17), and thereafter in T of the action S in Eqs.(19) 
and (25). This procedure is identical to how we treated the huid component in converting T 
of Eq.(18) to T of Eq.(19). Substituting into Eqs.(17), we have 

the new expressions for the Einstein equations: 

= l = + (34) 

where 14// = Vmf{q) + Pf + o); JPi'vf = Ng{q'Y/2. We also have a new expression for 

the ehective quintessence action: 

S,.e,, = 1^1 e^dT^2gV,„(l - Ng[q'fl2), (35) 

Thus, very diherent from what Eq.(19) or (25) may have suggested, the Lorentz factor here 
contains only the bulk flow, and the ’’thermal” energy contribution is absorbed into W//- 
The functional form of V depends on what kind of disk conhgurations under consideration. 
The new effect arising from will lead to a qualitatively very different result when the 
potential is attractive, where the nearly balance of the repulsive force of and the attractive 
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force of V gives arise to quasi-equilibrium trajectories, which on one hand oscillate on the 
fast time scale and on the other hand sink to the potential bottom on a secular time scale. 

It is of little surprise to hnd that the effective action of Eq.(35) also describes the dy¬ 
namics of a single complex held [^], the so-called ’’spintessence” where we let = 1 
and m = 2. The new effect pertaining to U^{q,a) is identical to the ’’centrifugal potential” 
of ’’rotational motion” of the single complex quintessence. This is because the ’’thermal 
motion” allows the trajectory of each individual held to carry a hnite ’’angular momentum” 
relative to other helds, and these transverse ’’velocities” give rise to the ’’centrifugal forces”. 
However, unlike ’’spintessence”, where the U{1) symmetry is globally destroyed to a large 
degree, the ’’thermal motion” of N helds break the U{1) symmetry to a much lesser de¬ 
gree. In particular, when the coarse-grained average of helds of a large N is taken, the U{1) 
symmetry is nearly restored, at worst broken by due to the Poisson statistics. 


D. Hamilton-Jacobi Equation 


We proceed to discuss how the present formulation can be extended to the Hamilton- 
Jacobi theory, which will prove to be useful for considerations of the cosmological pertur¬ 
bations of superhorizon size to be discussed later. We shall now return to the general 
expression, Eq.(18), for the many-held action. Given the Lagrangian in Eq.(18), we hnd the 
canonical momenta of the helds and huid to be 

m = gfij[^2Veff - Pf Tif = 9^{ - Pf -Pf)e^'^, (36) 

Pf+Pf 

where the Lagrangian multipier A has been set to its due value, 1. As the Hamiltonian 
h = it follows from Eq.(17) that 


= -Pf- Pf)e^^- 


n 


(37) 


/ 


We now employ the Hamilton-Jacobi theory to rewrite this Hamiltonian into 


-YW) -9{2V^ff-Pf-Pf)e , 


(38) 
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where the canonical momenta tt* and tt/ are cast into dS/d'ipi and dS/dOf respectively, 
the Hamiltonian h = gdS/dT and S is the Hamilton principal function. The so¬ 
lution to this Hamilton-Jacobi equation has the form S = Sqigpi^T) + Sf{9f), where 
Sq = Wq{'ipi,T) exp{3T). The fluid part Sf on the right-hand side of Eq.(38) cancels by 
itself since the time component of the fluid 4-velocity, nf9f/{pf +Pf), equals to unity, and 
it yields tt/ = y/guf exp{3T). Furthermore, the fluid momentum tt/ turns out to be a con¬ 
stant since the fluid solution obeys Uf oc exp(—3T), thus allowing Sf consistently to be 
independent of T explicitly. Equation (38) then becomes 

At any given time slice T = T„, the phase function Sq is a function of V’i, and dWq/dxpi 
is parallel to the velocity %!)[. Any virtual displacement lying on the 

Nm — 1 dimensional hypersurface Hdj(t/’i, T„) = const, is perpendicular to the local velocity 

where and are two adjacent instantaneous trajectories at T = Tq. At the next 
instant, T = Tq + 6T, the hypersurface evolves and so do the two trajectories and 
However, at this next instant the two trajectories may not land on the same constant- 
Wq hypersurface. In order for the two instantaneous orbits be on the same hypersurface, 
they must be evaluated at different times. The progressively asynchronous orbits on the 
constant-lTg hypersurfaces occur regardless whether the fluid is present. 

The notion of the constant-Hdj hypersurface is useful only when Aipi is sufficiently small, 
i.e., in the regime of linear perturbations. The reason is that as the multi-held dynamics 
is likely chaotic, Wq is generally not an analytical function of tjji. Generally speaking, the 
constant-lTq surface should consist of foliated patches jointed by kinks and cusps, and an 
analytical Wq exists only in some immediate neighborhoods of a given location tjji. Hence, as 
long as A'lpi is small, two trajectories can always stay in the same foliation patches at different 
times. In the next section, we will make use of this notion to discuss the cosmological linear 
perturbations in the comoving gauge. 
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IV. COSMOLOGICAL PERTURBATIONS 


The above formulation using T(= In a) as the time variable is a natural choice for describ¬ 
ing the coupled dynamics of the background scalar field and gravity. It turns out that such 
a choice is also convenient for analyzing the cosmological perturbations of space curvature 

0 . 

In the very long-wavelength limit, for which the wavelength k~^ is much greater than 
the Hubble radius H~^, the quintenssence and matter perturbations, and 6pf, can be 
regarded as spatially homogeneous. The energy-momentum tensor of each matter component 
is diagonal in the limit k/H 0, since the off-diagonal component involves two spatial 
derivatives and hence is of order {k/H^ and the space-time component is of order {k/H). 
Such a spatially homogeneous diagonal energy-momentum tensor leads to diagonal metrics. 
Thus, superposed with super-horizon perturbations of /c/iJ —> 0, both scalar held and huid 
satisfy the same equations of motion as those of the background derived in the last section. 
The cosmological perturbations are therefore the deviation of different trajectories of slightly 
different initial conditions. 

In fact, one can do better than the above zero-th order (of k/H) considerations and 
retains the perturbed quantities of order k/H. It is well known that the metric huctuations 
have gauge degrees of freedom, and some gauge choices may be more convenient than others. 


As has been shown by Sasaki and Tanaka [|I3, the gauge that keeps the time variable T 
unperturbed turns out to be the convenient gauge in allowing one to calculate perturbations 
of order k/H, specihcally for the perturbed space curvature R and the traceless stress. In 
this gauge, the perturbed shift (space-time component) of the metrics vanishes and the 
perturbed space compression (trace of space component) of metrics remains a constant. 
The equation of motion derived above, describing the ”00” component and the trace of 
the space components of the energy-momentum tensor, remains valid even up to 0{k/H) 
and is decoupled from other tensor components of size {k/H), namely, the perturbed space 
curvature and traceless stress. The great advantage of this gauge over, for example, the 
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synchronous gauge is that the perturbed space curvature can be directly determined from 
the freely decaying traceless stress. 

Although Sasaki and Tanaka’s work does not include the fluid component, the form of 
the metrics remains the same as the pure-scalar held case when the huid is present. One can 
show that the leading-order contribution from the huid in this constant T gauge is from the 
” 00” component and the trace of space components of the perturbed huid energy-momentum 
tensor. Moreover, as the perturbed Spf and Spf are related by the underlying equation of 
state, Spf is simply co{pf+pf) with a constant cq. In addition, the perturbed time-component 
of the huid 4-velocity SUo is simply —SH/H, the perturbed Hubble parameter; the space 
component SUi appears only on the order of 0{k/H) in the ”0i” component of the perturbed 
energy-momentum tensor in the same manner as the scalar helds do. Hence, inclusion of the 
huid retains all the good features of the constant-T gauge discussed by Sasaki and Tanaka 
for pure scalar helds. This is due partly to that the huid dynamics can be formulated in a 
way very similar to that of the scalar held, as has been shown in the last section, and partly 
to that the formulation using T as a new time variable has already absorbed the gravity in 
the matter dynamics in a self-consistent manner. 

Since the cosmological perturbations result from the departure of slightly diherent ini¬ 
tial conditions, there exist 2{Nm -|- 1) independent perturbations, due to the 2{Nm + 1) 
independent phase-space coordinates for the N scalar helds of m components, plus a huid 
component. After one manages to make the 2{Nm + 1) scalar modes orthogonal, the space 
curvature of the constant-T hypersurface only exists in two canonically conjugated modes, 
and vanishes in all others. These two modes are the adiabatic growing and decaying modes, 
and the rest are entropy modes. The growing adiabatic mode is associated with two trajec¬ 
tories with exactly the same initial phase-space coordinates but dihering by a ATq in their 
initial launching times. The two orbits follow exactly the same phase-space trajectories, 
but at any given time T their phase-space coordinates are diherent. This time-translational 
mode exists regardless of whether the huid is present. 

Among all independent perturbations, some may be integrable and thus stable, some 
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are chaotic and thus unstable, and the rest can still follow some stable manifolds in the 
phase space and thus be stable. Note that the time-translation mode is always stable. For 
those perturbations that are integrable, they are all associated with some integrals of motion 
ctj-’s, and result from different trajectories with slightly different a/s. For example, the above 
Coulomb field gas has a conserved total angular momentum, L, and there exists an integrable 
angular-momentum perturbation. Since the Coulomb gas repels, the ’’rotational motion” 
must slow down and hence the energy density of the angular-momentum mode decays in 
time. For those that follow stable orbital manifolds, these modes eventually converge into 
some attractors since this dissipative nonlinear system has attractors in various regions of 
the 2Nm dimensional phase space. These modes are usually associated with the ’’velocity” 
degree of freedom, as the presence of cosmic drag makes velocities eventually die out. By 
contrast, the chaotic unstable modes are usually associated with the ’’position” degree of 
freedom since the initial transient can render the orbits to end up at very different locations; 
these chaotic unstable modes give rise to perturbations in the potential energy V^ff- In great 
contrast to the single-held case, there are many independent entropy modes in the multi- 
held system, and many of these entropy modes can be more unstable than the adiabatic 
huctuations. 

As the multi-held dynamics is likely chaotic, the perturbations tend to be entangled. 
Though the constant-T slice provides a convenient way to calculate the cosmological per¬ 
turbations, the physically relevant time slice is such that the perturbed energy-momentum 
tensor has a vanishing Oi component, i.e., the comoving slice. This is due to that the 
constant-T gauge contains an extra gauge mode, which is always attached to the time- 
translational mode mentioned earlier. This gauge mode can be removed in the comoving 
gauge. 

The matter rest frame requires Tg = 0, or 

= 0, (40) 

* Pf + Pf 

where k is the small but hnite wavenumber of the space Fourier mode. It means that the 
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difference of two adjacent trajectories is constrained to be perpendicular to the background 
’’momentum” (c.f., Eq.(36)). As A'lpi and A6f lie on a Nm dimensional comoving hypersur¬ 
face in the Nm -|- 1-dimensional embedded space, this Nm dimensional comoving time-slice 
is simply the S{'ijji, Of)\T = const, hypersurface, where S is the Hamilton principal function. 

We recall the discussions at the end of last section that two adjacent orbits must be 
evaluated at different times AT for them to he on the same constant-Hd? hypersurface. 
Now, the arguments should be extended to include the fluid, and Wq is replaced by S. The 
difference of nearby trajectories 6ijj and 66f at the constant-S' surface therefore satisfy 

Y.i’HmT) + i>'AT) + + 0'jAT) = 0 , ( 41 ) 

i Pf + Pf 


thereby giving the time difference 


AT = 


+ {n}/pf +Pf)sefe'j 


(42) 


EMiy + {n}/pf + pf){9'fy 
where 6'ijj are 66f are the difference of two adjacent orbits evaluated at the same T. For the 
time translational mode, we have vanishing 6(j)i and 66f when observed in the comoving slice, 
and hence ATq = = 6df/d'p thereby yielding AT = —ATq. This cancels the original 

enhancement of curvature perturbation R due to the difference of initial launching times 
ATq in the constant-T gauge, and reverts the curvature perturbation back to its intrinsic 
value Rq given by initial quantum fluctuations. That is, the contribution of ATq to the 
perturbed space curvature R in this constant-T gauge is the gauge degree of freedom and 
removed in the comoving gauge 0,0 • 


V. CONCLUSIONS AND DISCUSSIONS 

We formulate the coupled dynamics of quintessence, gravity and fluid in a unihed frame¬ 
work by using the logarithmic cosmic expansion factor T(= In a) as the new time variable. 
The held dynamics turns out to be identical to the dynamics of a relativistic particle im¬ 
mersed in a static viscous medium. The presence of a huid component serves merely to 
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weaken the force acting on the held. Such formulation can be extended to describing many- 
held dynamics as a many-body problem. The forces on these helds can be made to be 
self-consistent ones by incorporating the mediate scalar bosons in the held functional space. 
This work also points out that when the number of helds is large, it is possible to reduce 
the many-held problem to an ehective single-held problem. We work out the cold Coulomb 
held gas and warm Coulomb held gas as two examples. 

It is interesting to note that when the multi-held dynamics is cast in the framework 
of relativistic Hamilton-Jacobi theory, it can be straightforwardly extended to quantum 
mechanics where Eq.(38) becomes the Klein-Gordon equation with an evolving squared 
mass, 2gVeffe^'^ . This can be made exact when pj = 0 in the very early universe before 
photons were produced. The quantization are with respect to the scalar held tjji relative 
to the conjugate momentum (c.f., Eq.(36)), as well as the Hamiltonian h relative to 
the conjugated expansion factor lna(= T) (c.f., Eq.(37)). Note that the gravity has been 
incorporated in the quantization since the action of Eq.(18) contains the gravity piece. Due 
to the presence a time-dependent squared mass in the Klein-Gordon equation, it gives rise 
to low-energy (~ hH) quantum huctuations of the quintessence held and metrics on the 
superhorizon scale. 

The quantization of tjji in the held space logically hts well with our earlier treatment 
of the scalar bosons V, which mediate the quintessence interactions discussed in Sec. (3), 
in the held functional space. (Of course, quantization of the scalar bosons V can also be 
performed when necessary.) As has been pointed out immediately below Eqs.(7) and (8), the 
ultra-relativistic regime renders the Hamiltonian h to be a good constant of motion. In this 
regime, the interactions Vij is comparatively small, and the dominant energy component is 
the free held. It thus provides an easy way for quantization of ipi together with the gravity 
along the light cone. 

The space coordinates have been dropped out of scene in the above analyses because we 
are conhned to the superhorizon scales. It may be regarded as a particular reference frame 
where the dynamics takes place in the rest frame of the space coordinates (in contrast to the 
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field coordinates) with a vanishing space momentum. It remains to be investigated whether 
the Planck-scale physics really admits such a symmetry that one can generally transform 
the reference frame to the proper frame where the space momentum vanishes. For example, 
when the uncertainty principle applies to the classical notions of space and time in the Planck 
era, space can become fuzzy in a reference frame where time has a definite value. If so, the 
Bose-Einstein condensation of a collection of N such superhorion quanta may constitute the 
foundation for the background N quintessence fields as well as the gravity in the Planck 
scale, and provide the needed energy content to develop into one universe. 

We now return to discussions on the classical quintessence helds. Clearly, the dynamics of 
classical multi-field can easily be non-integrable and chaotic. Though a large number of helds 
help simplify the problem, such a step is probably too drastic to take for a possible extension 
beyond the single-held quintessence. It is probably more natural to investigate dynamics of a 
few coupled helds together with a huid component. However, this is surely a difficult regime 
as far as the chaotic dynamics is concerned, since the coarse-grained averaged quantities 
do not exist. Nevertheless, since the friction applies uniformly in the velocity space, the 
acceleration may quickly become negligible and the velocities become non-relativistic after 
an initial transient, especially when the huid energy density is non-negligible and weakens 
the forces on the quintessence helds. In such a ’’rolling” regime, the degrees of freedom are 
reduced by one half and the dynamics becomes much simpler. However, even after such 
reduction, if the remaining degrees of freedom exceed or equal 3, it is possible that there 
exists strange attractors, where the trajectories stretch, twist and fold in a multi-dimensional 
manifold, and one will again face a technically insurmountable problem. This can be so when 
the number of coupled quintessence helds exceeds or equals 3. 

On the other hand, the scalar-held dynamics in the early universe has gross imprints 
in the late-time cosmological perturbations. Some entropy perturbations may have larger 
held amplitudes than the adiabatic perturbations do, if indeed multi-helds were active in the 
early universe. They generally yield a diherent huctuation power spectrum in the cosmic 
microwave background radiation (CMBR) from the one predicted by the pure adiabatic 
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perturbation. To be consistent with the observed CMB temperature anisotropy spectrum, 
/c” with n K, 1 over three orders of magnitude in length scale, or in unit of the logarithmic 
scaling factor AT ~ 10, the entropy modes cannot be too unstable within this interval in 
the inflationary epoch. This is due, on one hand, to that the potential V is flat during 
inflation and the field dynamics is slow, and on the other hand to that the entropy modes 
normally lead to a spectrum of density perturbations dominated by by the large scale, but 
the COBE measurements have fixed their amplitudes to be less than that of the adiabatic 
perturbation |^. However, the entropy mode may become unstable after the inflation. 
Instabilities give rise to small-scale power in the CMB fluctuations. In view of the relatively 
large error bars in the CMB measurements of fluctuation spectrum [^|-|^ , there may still be 
rooms for accommodating entropy modes resulting from the multi-field quintessence. 
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